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Abstract 

The spontaneous emergence of enhanced responses and local orders are properties often associ- 
ated with complex matter where nonlinearities and spatial inhomogeneities dominate. We discuss 
these phenomena in quantum devices realized with superconducting Josephson junction networks 
and cold atoms in optical lattices. We evidence how the pertinent engineering of the network's 
shape induces the enhancement of the zero-voltage Josephson critical currents in superconducting 
arrays as well as the emergence of spatially localized condensates for cold atoms in inhomogeneous 
optical lattices. 
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I. INTRODUCTION 



In all fields of physics homogeneous systems have the simplest properties and, thus, play 
a very important paradigmatic role in our understanding of natural phenomena. It is a fact, 
however, that many real systems are inhomogeneous in a way or another and that, in some 
instances, their inhomogeneity may be the seed for the emergence of new and unexpected 
complex behaviors [1], which may be probed in experiments and, hopefully, become useful 
in the engineering of quantum devices. 

In condensed matter systems, inhomogeneities may lead to enhanced responses to external 
perturbations and/or to the emergence of local orders 2|]. Remarkable examples include the 
large transport anisotropy observed at low-temperature in quantum Hall samples after the 
onset of electronic nematic phases [31], the colossal magneto- resistance in manganites J], 
the appearance of striped phases in systems with competing interactions as well as of 
pseudogap phases in high superconductors ^. In all the above mentioned examples 
the relevant optimal inhomogeneities are dynamically generated l7|]; furthermore, stripe and 
pseudogap phases are associated to the onset of local orders js], [J. 

To get control on the onset of complexity of a condensed matter system is desirable 
not only for understanding new emergent functionalities useful for the engineering of new 
materials and devices but also for discovering levels of theoretical description enabling to 
well separate the properties of global phases from the ones arising from phase competition; 
global average behaviors are, in fact, not helpful for this task |9|]. Quantum devices pro- 
vide a controllable setting (one can fabricate them!) to investigate the effects induced by 
inhomogeneities on the emergence of complex behaviors. 

Quantum devices with built-in inhomogeneities may be realized with today's available 
technologies using either superconducting Josephson junction networks (JJN) lOj or ultra- 
cold atoms in optical lattices pJL|]. Due to their versatility and to the great reliability of 
the fabrication technologies developed for their construction, JJNs and ultracold atoms in 
optical lattices are by now the prototype of complex physical systems exhibiting a variety 
of interesting physical behaviors, adjustable by acting only on very few external parameters 
and, as we shall point out, also by the pertinent engineering of the network's shape; in 
addition, they provide controllable settings to investigate the properties of granular super- 
conductors or high Tc superconductors [12] paving a very promising avenue in the engineering 
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of quantum states of potential interest for processing quantum information 13|. 

A Josephson junction may be realized using two superconducting grains separated by 
an insulating layer; by an appropriate engineering of the insulating support, one can fabri- 
cate an inhomogeneous network of Josephson junctions with a given shape. Inhomogeneous 
JJNs [isl have been studied for a long time with the aim of establishing the effect 
of space connectivity on superconductivity |l6]. Recently, the appealing perspective to 
realize devices for the manipulation of quantum information stimulated the analysis of in- 
homogeneous planar JJNs with non conventional connectivity 17(1, engineered to sustain a 
topologically ordered ground state Furthermore, transport measurements on supercon- 
ducting wire networks evidenced - even in pure systems with non-dispersive eigenstates - 
interesting anomalies of the network critical current induced by the interplay between the 
network's geometry and topology with externally applied magnetic fields [1.9^. Furthermore, 
the theoretical analysis of rhombi chains has evidenced the exciting possibility of being able 
to detect 4e superconductivity thr oug h measurements of the supercurrent in presence of a 
pertinent external magnetic field 201]. Here we shall address the properties of inhomoge- 
neous JJNs fabricated on lattices with non-random regular inhomogeneities engineered to 
yield enhanced zero-voltage Josephson critical currents as well as local order on domains 
selected by the network's topology. We shall investigate in detail the paradigmatic example 
of a comb-shaped JJN, whose properties have been analyzed in 2l|, |22]. To fix the ideas, 
in Fig. [T] is reported the design of the device; there, the circles locate the position of the 
superconducting grains while the crosses represent the junction joining them. 

Control over the network connectivity may be achieved also with ultracold atoms on opti- 
cal lattices [23|. A one-dimensional optical lattice is realized using two counter-propagating 



laser beams arranged to create a periodic potential where bosons may tunnel between wells 
with a rate adjustable by tuning the laser's power controlling the height of the inter-well 
barrier. In this realization a lattice site is a minimum of the periodic potential while the links 
between neighbor vertices are provided by the barrier between wells. By suitably arranging 
more laser beams one is able to engineer the shape of the optical networks allowing, in 
principle, for an experimental testing of the proprieties of inhomogeneous lattices such as 
the ones analyzed in 
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26|. To fabricate bosonic Josephson junctions (BJJ) and networks 



(BJJN) one should require also that all the atoms in a given well are described by the same 
macroscopic wave function: a BJJ is obtained, then, by loading - at a temperature below 
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FIG. 1: The comb graph. 



the Bose-Einstein condensation critical temperature - a condensate in each well of a double 



we. 



1 potential [27[ while the barrier separating the two condensates acts as a Josephson link 
281] . With a multi-well optical potential, provided that the heights of the barriers are much 
higher than the condensate's chemical potential, one fabricates a BJJN (i.e., a lattice of 
weakly coupled condensates). In the following we shall evidence how bosons hopping on 
comb-shaped optical lattices undergo a spatial Bose-Einstein condensation on the comb's 
backbone; this happens even if the bosons are free and the network' s euclidean dimension 
is 1. 

The emergence of complex behaviors in the superconducting and bosonic arrays analyzed 
in this paper is determined solely by the pertinent engineering of the network's connec- 
tivity, since it is due only to the spectral properties of the adjacency matrix (i.e., of the 
matrix defined as Aij = 1 if i and j are vertices connected by a link, and Aij = otherwise 



291]), which fully characterizes the network's geometry and topology. In some very specific 



instances 
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3l| , the spectrum of Aij contains a continuous set of states - the hidden spec- 
trum - localized around domains selected by the network's topology with eigenvalues ranging 
continuously from a certain value Eq up to the threshold of the continuum delocalized states 
Efi'. the spectrum is then effectively gapless even if, in the thermodynamic limit, there is a 
lowest eigenvalue Eq, confined away from E^. In the following we shall evidence that it is this 
spectral anomaly - induced from the pertinent choice of the network's connectivity- which is 
responsible for the emergence of enhanced responses and local orderings in superconducting 
and bosonic networks, respectively. 



An hidden spectrum of the adjacency matrix emerges, for instance, when one analyzes 
bundled graphs {29! (i.e., those obtained by grafting a fiber graph to every point of a base 
graph) while, for graphs with constant coordination number (such as the Sierpinski ^ gasket 
and the ladder graph), the adjacency matrix does not support any hidden spectrum [31]. In 
the following, we shall analyze the simple paradigmatic case of comb networks showing how 
the hidden spectrum of the adjacency matrix leads to unusual quantum behaviors such as 
the emergence of the spatial BEC on the comb's backbone for a Bose gas living on a comb- 



shaped optical lattice [30|, l32| and of the enhanced responses observed for classical comb 



shaped JJNs made of Niobium grains |2lL |22[[see Fig.l]. To better clarify our arguments. 



we find instructive to compare our results with those obtainable if the same devices were 
defined on a chain, since the latter is, after all, the simplest graph of euclidean dimension 1. 

In section 2 we determine the equations governing the properties of inhomogeneous su- 
perconducting and bosonic networks. Section 3 is devoted to the analysis of the spectrum 
of the adjacency matrix of a comb-shaped network and to a comparison of its spectrum 
with the one pertinent to a linear chain: there we see that, while the adjacency matrices of 
both graphs admit a continuum set of delocalized states starting at the same eigenvalue, the 
adjacency matrix of a comb supports an hidden spectrum. In section 4 we evidence how the 
states belonging to the hidden spectrum are responsible for the enhancement of the zero- 
voltage Josephson critical currents in superconducting JJNs as well as for the emergence 
of the spatial BEC in bosonic networks on domains pertinently selected by the network's 
topology. Section 5 is devoted to our conclusions and final remarks while the appendices 
report on some pertinent but rather lengthy computational details. 



II. MANY-BODY THEORY OF SUPERCONDUCTING AND BOSONIC SYS- 
TEMS ON GENERIC NETWORKS 

In this section we shall summarize the many-body description of superconducting and 
bosonic networks. Our purpose here is mainly to clarify the assumptions underlying the 
equations used in section 4 to characterize the complex behaviors emerging from a pertinent 
choice of the connectivity in JJNs and bosonic networks. Our approach uses the self- consis- 
tent Bogoliubov-De Gennes (BDG) equations ^] since they provide an unified framework 
to account for the description of fermionic and bosonic systems enabling to appreciate in a 
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rather simple context the differences (and similarities) between these systems. 

We shall describe first the microscopic theory of inhomogeneous superconducting net- 
works realized by putting on each site of a graph a superconducting grain. On each island, 
the effects of the electron-phonon and Coulomb interactions are embodied in the BCS pa- 
rameter and, furthermore, there exist a critical temperature Tbcs such that for T lesser 
than Tbcs each grain becomes superconducting; of course, the entire network becomes su- 



perconducting at a lower temperature Tc 



34| . We, then, analyze ultracold atoms (spinless. 



for simplicity) in deep optical lattices and derive the set of equations describing both the 
excitations and the condensate's dynamics; in dealing with atomic systems one should only 
account for repulsive interactions between bosons as it is well known that attractive inter- 
actions lead to instabilities even in presence of confining traps 23 1. 



A. Superconducting networks 

Inhomogeneous fermionic systems with attractive interactions may be conveniently de- 



scribed by the BdG equations 



331]: with a two-body point- like interaction V{f — f) 



-V(5(r — f), the Hamiltonian is 7-^ = Ho + Hi, where 

Ho = j dr^m\fa)Ho^{fa) (1) 



and 



7^1 = -^ j df^¥{ra)-^\ra')-^{ra')-<^{fa). (2) 



In Eqs. (11]) and ([2]), the \l/'s are fermionic operators, cr = ± is a spin index, and Hq = 
— /i^V^/2m + [/o(r') — fi, where /i is the chemical potential and f/o(r) is an external potential 
(we assume, as usual, that f/o(r) is spin-independent and that no magnetic field is applied). 
When the fermions are not interacting, then fi = Ep, where Ep is the Fermi energy. In a 
self-consistent approach, one defines the effective Hamiltonian 



ra 



+ y"c/f|A(f)^"f(f+)^t(f_) + A*(f)^(f-)^(f+)| (3) 

and, then, requires that the pair potential A(r) = V(\l'(r+)\E'(r— )) and the Hartree-Fock 
potential U{f) = — V(\l'^(rcr)\I^(fcr)) are self-consistently determined from the solutions of 



the BdG equations 

G„ Ua{r) = [Ho + f/(r)]u,(f) + A(f)i;,(f) (4) 
e„ v^ir) = -[Ho + Uif)]v^ir) + A*(f)n„(r), (5) 

with Ua, Va satisfying 

J df[\u^{f)\' + \v^{f)\'] =1. (6) 

One has then 

A{r} = v5^M,(rl<(r)tanh e^) (7) 

and 

U{f) = -Vj2 [Mr-)\'fa + Mr-)\' (1 - /.)] , (8) 

a 

where /„ = (e^^° + (3 = l/ksT and the sums are only taken on positive Gq,'s. The 
chemical potential fi - when different from the Fermi energy - is determined from the normal- 
ization condition. Of course, when there is an external potential breaking the translational 
invariance, the pair potential is dependent on the position. 

For spatially homogeneous superconducting networks the quantum number a is just the 
momentum k and, then, A = A{r) and U = U{f) do not depend on the position implying 



that Gg= y A^ + E|, where = h'^k'^ /2m — jj + U . Furthermore, since u^{r} = L ^/"^U^e^^'^ 
and v^{r) = L^'^/'^V^e^^"^ with being the volume, one gets U'i = (1/2) [1 + Ej^/ G^] and 
= (1/2)[1 - E^/ Gg]. With the BCS interaction (i.e., V = Vbcs if \h^ky2m - Ef\, 
\h'^k''^/2m — Ef\ < hujo and V = otherwise, ujd being the Debye frequency) and f/ = 
(which is a good solution of the self-consistent equation for U if n(0)VBcs ^ 1), one has 
H = Ep together with the celebrated BCS gap-equation 

2 VA' + E^- \2 }' ^' 

n(0) is the density of states per volume and spin direction at the Fermi energy. For inhomo- 
geneous networks, instead, one should regard the external potential Uo{r) as representing 
the insulating barriers between the grains and take self-consistently into account its effects. 
To obtain a discrete version of the BdG equations (LBdG) suitable to describe the classical 



superconducting JJNs fabricated in 21|, one makes the ansatz that the eigenf unctions of 



the continuous BdG equations 331] may be written in a tight binding form as Uair) 



^i^a(O0i(r) Va{r) = J2i'^a{'i)4>i{'f^] ^ labcls the position of a superconducting island 



while the contribution of the electronic states participating to superconductivity in a given 
island is effectively described by a field whose specific form depends only on the 

geometry of the islands and on the fabrication parameters of the connecting junctions. The 
assumption that 0j(r*) does not depend on a amounts to account only for contributions 
coming from electrons near the Fermi surface. The LBdG equations then read 

Ga Ua{%) = ^ tijUaii) + A{l)Va{i) (10) 
j 

j 

where Ua and Va satisfy to 

+ =1- (12) 

i 

The matrix e.^ is defined by 

eij = -tAij + U {i)6ij - fi5ij, (13) 
with Aij being the adjacency matrix of the network, 

fi = li- j df(f)i{r} {-h^V'^/2m) (f)i{r} (14) 

and 

t = - j df4>i{r)[-h^Vy2m + Uo{r)]4>j{r). (15) 

Self-consistency requires 

A(z) = V^n,(z)<(^)tanh (16) 

and 

uit) = -vJ2 [Km' fa + Km' (1 - /.)] , (17) 

a 

where V = V0^(r* = fi) is assumed to be independent of i. The network's topology and 
geometry is encoded in the term —tAij appearing in the definition of the matrix e^j given in 
Eq. f|T3l) . while the specific values of t and V depend - as a result of our ansatz on the form 
of the eigenfunctions of the BdG equations- only on the 0i(r). 

To justify the assumptions involved in the derivation of Eqs. ( fTOjl and (fTTjl . one should 



observe that, for the JJN device described in [2l|, capacitive (inter islands and with a ground) 
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effects are negligible, that the total number of electrons on the island TV is much larger than 
the number of electrons tunneling through the Josephson junction and that all the islands 
contain approximately the same M {Af{i) = A/"). Furthermore, the islands are big enough 
to support the same superconducting gap of the bulk material; as a result, one may require 
0j(r) to be position- independent on each island except for a small region near the junction 
and to be the same on each island with a normalization given by / df(f)'j{r) = N'{i) = M and 
j df(f)i{r)(f)j{r) ^ for i ^ j; in our derivations we set A/" = 1. As a result t ^ Ej = {h/2e)Ic, 
where Ej is the nominal value of the Josephson energy of all the junctions in the network 
while Ic is the bare zero-voltage Josephson critical current of each junction. In section 
IV we shall provide an explicit solution of the BdG equations (IIUI) and fill I) describing 
superconducting JJNs fabricated on a comb-shaped insulating support [21], while in App. 
A a solution of the LBdG equations for a Josephson junction chain is provided. 



B. Bosonic networks 



The full quantum Hamiltonian pertaining to a (spinless) bosonic gas in an optical poten- 
tial is H = Ho + "^i, where 



(18) 



and 



Hi 



go 

2 



df<^\r)<^\r}<^{r}<^{r). (19) 

In Eqs. ( |T8l) and ( |T9l) . $ is bosonic operator, and Hq = — ^^V^/2m -|- Uo{r) — fi; Uo{r) is the 
external potential, which is usually the sum of the optical lattice and of the magnetic trap 
potential; for simplicity, we do not consider terms accounting for the effects of harmonic 
traps. In writing Eqs. ( fTSl) and ( fT9|) . the standard s-wave scattering approximation has 
been used: i.e., the two-body potential is written as l^(r — r') = go6{f — r'), where go = 
Anh^a/m > 0, with a being the s-wave scattering length and m the atomic mass. The 



dynamics of the bosonic field is described by the well known equation 



23\ 



iTi—^ 
dt 



(20) 



from which the Gross-Pitaevskii equation for the condensate wavefunction = ($) is usually 



derived 



23|. 
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The bosonic counterpart of the BdG Eqs. (jlj) and ([5]) is provided by the Bogohubov 



231 ]: to derive them, one usually writes $ = tp+S^ 



equations for the condensate's excitations 
with 

5$ = 5^K(F)a„e-^^'^*/'* - vl{r}aie'^-'/% (21) 

a 

where are operators destroying bosons in the excited state a. Using Eq. (!20|) and keeping 
only terms linear in the fluctuation field 5$, one finds 

Ea Ua{r) = [Ho + 2gono{r)]ua{r) - gono{r)va{r) (22) 

Ea Uair) = [-Ho - 2gonoir)]uc,ir) + gono{r)Vair) (23) 



where no{r) = |'?/'(^)P is the condensate density [23|. Notice that Ua and satisfy now to 
the condition 

j df[\u^{f)\'-\v^{r)\'] = l. (24) 

In the analysis of a bosonic system one has to use, in addition to the BDG equations, 
the Gross-Pitaevskii equation for the condensate. To do this, one should observe that, when 
the power laser is high enou gh, the bosonic field $ may be approximated by means of the 
tight-binding approximation [35| as 

mt) = Y.b,mi^-^ (25) 

substituting this ansatz in the full Hamiltonian, one gets (in the non interacting limit) a 
simple tight-binding model described by 

H = -tJ^A,,blb,, (26) 

where the coefficient t is given by Eq. ( ITSil . In Eq. ( l26il . i,j denote the minima of the optical 
lattice (i.e., sites of the network) and 6^ (bj) is the bosonic operator which creates (destroys) 
a boson at site j. The filling, i.e., the average number of particles per site, is defined as 
/ = Nt/Ns, where is the total number of bosons and Ns is the total number of sites. 

Eq. (|26l) is the pertinent equation to investigate in order to ascertain if , for a gas of 
ultracold bosons, BEG emerges as a result of the network's inhomogeneity. When BEG 
occurs in each well, each pair of neighbouring wells acts as a bosonic junction with Josephson 
energy given by 

Ej ^ 2tf. (27) 
10 



In the next sections we shall show that a condensate indeed emerges when bosons hop 
on comb- shaped networks. 



III. SPECTRUM OF A QUANTUM PARTICLE HOPPING ON COMB-SHAPED 
NETWORKS 

In section 2 we evidenced how the equations describing superconducting and bosonic 
networks depend on the adjacency matrix characterizing the network's connectivity. In this 
section we shall review [s^ the main results concerning the spectrum of the adjacency matrix 
describing the connectivity of a comb network. 

A comb (see Fig. [1]) is made of one-dimensional chains {fingers) grafted periodically on 
a linear chain (backbone). Each site of the comb can be naturally labeled by introducing 
two integer indices (x, y), where x = 0, - ■ ■ ,Li labels the different fingers and y = 0, - ■ ■ , L2 
provides the distance from the backbone. Each site on the finger is linked to two neighbors 
whereas each site of the backbone has four neighbors. 

The topology of the network is fully described by the adjacency matrix A^^y.^'^yi which 
equals 1 if {x,y; x',y') is a link and otherwise. A quantum particle hopping on a comb is, 
then, described by the Hamiltonian 

H = —t Ax^y-x',y' bl ybx^yi- (28) 



The single-particle energy spectrum is found by solving the eigenvalue equation j30|, l3l|, l32| : 

- t ^ A^^y. ^/y ^Jje{x', y') = EipEix, y); (29) 

x' ,y' 

on a comb Ax^y-x',y' is given by 

Ax,y;x',y' = {Sx,x'+1 ,x'-l)Sy,0^0,y' + i^y,y'+l + Sy,y' -l)Sx,x' ■ (30) 

In the following, we shall determine the spectrum of a quantum particle on the finite Li x L2 
comb and, only at the end, take the limit Li,L2 00. On a finite Li x L2 comb - using 
the adjacency matrix ( l30l) - the eigenvalue equation ( l29l) reads: 
Li-l L2-1 

- ^ ^ ^ [{Sx,x'+1 + Sx,x'-l)Sy,oSo,y' + {Sy,y'+1 + Sy,y'-l)Sx,x'] , y') = El/j^X, y). (31) 

x'=0 y'=0 

11 



Without loss of generality one may take x and y to be positive integers, since, due to periodic 
boundary conditions, (0, 0) = (Li, 0) and (a;, 0) = (x, L2). The total number of sites is then 
Ns = Lix L2. 

By exploiting the translation invariance along the backbone, a Fourier transform in the 
variable x reduces Eq. f l31l) to a one- dimensional eigenvalue problem for a quantum particle 
hopping on the comb's fingers. In fact, upon defining 

^(A;,y) = 5^e^'=^^(a;,y), (32) 

X 

with k = 2TTn/Li, n = 0, . . . Li — 1, the eigenvalue equation (1311) becomes: 

Li-l L2-I 

- t ^ ^ [2 cos{k)6yfl6o^y>Sk,k' + {Sy,y'+i + Sy,y'-i)h,k'] ^(^', v') = Etp{k, y). (33) 

k'=0 y'=0 

Since Eq. (1331) is diagonal in k, it may be written as: 

-t^[2 cos{ko)6y^oSo,y' + (^s/y+i + ^y,y'^i)] i'i.y') = E^piy), (34) 
y' 



where ilj{k,y) = 5{k — kQ)il){y) with fco = 27m/Li, n = 0, ...,Z/i — 1. Equation 
describes then a one-dimensional quantum particle interacting with a potential located on 
the backbone, Viko) = — 2tcos(A;o). 

To determine the eigenvalues and eigenvectors of Eq. fl34l) one may look for solutions of 
the form 

iPly) = Ao cos{hy + a) for -2t < E = -2t cos(/i) < 2t] (35) 
ipiy) = A_ e-^y + 5_ e'^y for E = -t(e'* + e"^) < -2t; (36) 
^{y) = A+{-iye~^y + B+i-iye^y for E = t{e^ + e'^'') > 2t. (37) 

To fix both the free parameters and the eigenvalues E one requires ip{y) to be normalizable 
and to be a solution of the eigenvalue equation in y = and y = L2 — I. These points are 
the only ones where Eq. (!34l) is not identically satisfied, yielding two equations to determine 
the two free parameters. Since for a given value of ko there are L2 different eigenvalues, the 



spectrum will consist of Li ■ L2 states and it can be divided in three regions: cxo and a± |30|. 
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(To is the part of the spectrum corresponding to delocahzed states with energies between —2t 
and 2t. Requiring the wavefunction fl35|) to be a solution of the eigenvalue equation flMl) in 
y = L2 — I and y = yields 

cos[h{L2 — 2) + a] + cos(a) = 2 cos(/i) cos[/i(L2 — 1) + a] 

(38) 

cos[h{L2 — 1) + a] + cos{h + a) + 2 cos(A;o) cos(a) = 2 cos{h) cos(q;), 
implying that there are odd eigenfunctions with a = 7r/2, h = 2Txm/L2 and m = 
1, . . . , L2/2 — 1 and L2/2 even solutions, obtained by 

- cos(fco) cot(/iL2/2) = sin(/i). (39) 

Equation (1391) can be solved graphically. In the large L2 limit the allowed values for h are: 
h ^ 7r(2m — 1)/L2, with m = 1, . . . , L2/2. For each value of ko there are L2 — 1 eigenvalues 
of type fl35|) with energy = —2tcos{h) and wavefunctions ilj{x,y) = e^''°^ sin{hy) and 
ifj^x^y) = e^^°^ cos{h\y \ + a), with /i and a satisfying Eq. (138|) . Thus, the fraction of states 
in this spectral region is / = Li{L2 — ^)/Li ■ L2. Of course, / tends to 1 in the limit Li, 
L2 ^ 00 and the density of states is given by |3l| : 

dE 

Po{E) dE = dn = L,d{L2h/27i) = Li{L2 - 1) , , ^ . (40) 
just as for a particle hopping on a linear chain. 
• cr_ 

(j_ is the part of the spectrum corresponding to localized states with energies E < —2t. 
Requiring that the wavefunction fl5Bl) is a solution of the eigenvalue equation flMl) in y = 
L2 — 1 and y = yields now 

A_ e-^(^^-2) + 5_ e^^^'-^^ + (A_ + fi_) = (e^ + e-^)(A_ g-'^^^^-i) + B_ e'^^^^"^)) 

^ ^ (^^e-/' + 5_e?») + 2cos(A;o)(^- +5_) = (e^ + e-'*)(A_ + fi_), 

(41) 

leading to: 

cos(A;o) coth(/iL2/2) = sinh(/i). (42) 

Equation fH2|) can be solved graphically, yielding a real solution only if cos(fco) > 0. The 
density of states is then given by 

\E\dE 

p_ (E) dE = dn = d(Liko/2n) = Li ' ' ^ . (43) 
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(T+ is the part of the spectrum corresponding to locahzed states with energies E > 2t. The 
parameters and -B+ are fixed by requiring again that the wavefunction fl37l) is a solution 
of the eigenvalue equation flMl) in y = L2 — 1 and y = 0. This requirement yields a set of 
equations similar to those of Eq. fj42|) . leading to cos(A;o) coth(/;,L2/2) = — sinh(/;,), which 
supports a real solution only if cos(/co) < while, for cos{ko) = 0, one gets the constant 
solution ip^y) = 1 with energy E = 0. The density of states in this spectral region is given 
by 

\E\dE 

p+(E)dE = dn = d(Liko/27i) = Li ' ' ^ . (44) 

The union of o"_ (i.e., states with E < —2t) and cr+ {E > 2t) forms the hidden spectrum 
and is the part of the spectrum corresponding to localized states; each spectral region 
contains Li/2 states and, thus, Li states belong to the hidden spectrum. Taking the limit 
L2 and Li 00, from Eq. ( l42l) . one easily verifies that the states of the hidden spectrum 



satisfy sinh(/i) = cos(/co) and correspond to energy eigenvalues E = —2t^\ + cos^ikt)) for 
cr_ and E = 2t^\ + cos^iki^ + 7r/2) for cr+. As the comb's size gets bigger, almost all the 
states- i.e., all the states apart from a set of measure zero - belong to do; in fact, since 
J^^Po{E)dE = Li ■ (L2 — 1), the normalized density of states belonging to ctq tends to 1 
while, for what concerns the hidden spectrum, one has lim2,i,L2-»oo /^^ P±{E) dE = 0. Since 
the hidden spectrum does not contribute to the normalized density of states p{E) of the 
pure hopping model on the comb, one has: 

Normalizing the density of states of the lower hidden spectrum to Li, one obtains for E E cf^: 

^ P- (E) = , = (46) 

An analogous equation holds for the spectral region E E a^. The density of states can then 
be plotted as in Fig. [21 where the pertinent normalizations for the continuous and hidden 
part of the spectrum have been used. 

The lowest energy eigenvalue in the normalized density of states is Ed = —2t] but, this 
is not the lowest energy attainable by a particle hopping on a comb since there is a lowest 
localized eigenstate belonging to cr_ whose eigenvakue is given by: 

Eq = -2V2t. (47) 
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FIG. 2: The density of states of the Hamiltonian (j28|) in units of t. The sohd line indicates the 
continuous part of the spectrum normahzed to Li • (L2 — !)• The dashed Unes denote p±, 
normahzed to Li. 



In fact, for L2 — ^ 00, since the energy is a decreasing function of h, the lowest energy level 
of cr_ is attained when cos(fco) = 1- 

Of course, Eq < Ed and this should indicate that the spectrum is gapped. However, one 
does not find an energy gap between Eq and Ed, since, for each value of (cos(A;o) > 0) 
there is a solution of Eq. (H2l) with a different energy in the interval [EQ,Ed\. In a finite 
comb of Li X L2 sites there are //i/2 solution of this type and, as Li 00, these solutions 
fill densely the interval [EQ,Ed\. These spectral proprieties should be contrasted with the 
ones arising when non interacting quantum particles are trapped in a harmonic well. 

From Eqs. ( l36l) and ( 14T1) it is also possible to show that, when Li, L2 00, the eigenvector 
corresponding to the lowest energy eigenvalue is 

i^E,{x,y) = ^,e-\y\'^, (48) 

where ^ = 1/h = 1/ log(l + \/2) is the parameter accounting for the localization around the 
backbone. In Fig. Owe plot the ground-state wavefunction as a function of the distance from 
the backbone; the plot well evidences the exponential localization only around the comb's 
backbone. 

Although not explicitly imposed, the condition of the continuity of the discrete gradient 
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FIG. 3: The normalized ground-state wavefunction of the comb network as a function of the 
distance y from the origin for a comb with 51 x 51 sites. 



holds at the points (0, y). Of course one cannot ask for the continuity of the linear derivatives 
along X or y since it is just this discontinuity which allows a particle moving on the finger 
to hop in the direction of the backbone. 

A similar analysis may be carried to determine the spectrum of a quantum particle 
hopping on different comb-like networks, such as the star-comb, mini-comb, and semi-comb 
depicted in Fig. |H It is not difficult to convince one-self that, for a pertinent choice of 
boundary conditions, the adjacency matrix describing the connectivity of these bundled 
graphs admits also an hidden spectrum. 

In the following we shall report only the lowest ener gy e igenvalues pertaining to quantum 



particles hopping on these networks. For a star-comb 36| one has 



Ei'=-tSlzf^, (49) 

where p is the number of arms on each star. For the semi-comb, one finds 

Eo^^™ = -2v^t, (50) 
where = (1 + v^)/2 ~ 1.618 is the golden section. For a mini-comb, one finds 



T-im 
^0 



- (l + v^TT^) t. (51) 
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FIG. 4: Top: The star-comb network. Bottom Left: Two "mini"-comb networks (top: p = 1; 
bottom: p = 2). Bottom Right: the "semi"-comb network. 



As we shall see in the following, the value of the ground-state energy is intimately tied 
to the critical temperature at which quantum bosons hopping on comb-networks undergo a 
spatial Bose-Einstein condensation. 

IV. COMPLEX BEHAVIORS EMERGING FROM THE NETWORK'S CONNEC- 
TIVITY 

In this section we review how the pertinent choice of the network's connectivity leads to 
;he emergence of new phenomena in quantum devices realized with superconducting JNNs 



22] and cold atoms in optical lattices 



30 



32[ |. Our subsequent analysis well evidences that 



new emerging phenomena are possible only if the network's connectivity is described by an 
adjacency matrix supporting an hidden spectrum. As we shall see, the states belonging to 
the hidden spectrum induce in a discrete many body system the onset of a new relevant 
energy scale and enhance the number of states which can be occupied by quantum particles 
at low energy. 
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A. Enhanced Josephson critical currents in a comb-shaped JJN 



In the following, we use the LBdG equations derived in Section 2 to compute the zero- 
voltage Josephson critical currents of junctions located on Josephson linear chains and 
comb-shaped Josephson networks. Using the eigenfunctions of the LBdG equations, a self- 
consistent computation yields for both systems the gap function, the chemical potential 
and the quasi particle spectrum. Our analysis evidences that, on the backbone of a comb- 
shaped JJN, the BCS equations are satisfied with a renormalized value of the Josephson 
energy. Then, we compute the zero-voltage Josephson critical currents Jc on the comb's 
backbone and compare our results for I^. with the outcomes of the experimental measure- 



ments presented in 



21 



22| and summarized in Fig. 5. 
For a comb network with L x L islands (see Fig.l), one finds a solution of the LBdG 
equations (fTOj) and (fTT]) where both the Hartree-Fock potential U{i) and the gap function 
A{i) are position dependent. The eigenvalue equation 

- AijlpaiJ) = Calpaii), (52) 

j 

leads to the emergence of the hidden spectrum analyzed in section 3. 

For a crude analytical estimate, one may require that, away from the backbone, the fingers 
may be regarded as a linear chain with uniform potentials (i.e., A{i) = and U{i) = Uc)- 
To get coupled equations for Af,, Ac, Ub, and Uc, one writes the LBdG equations (fTOj) and 
f fTT]) on a generic backbone's grain i. Upon setting Ua{i) = Uatpaii) and Vaij) = Va'ipaii)-! 
with f/^ + = 1, the self-consistency equation for U implies that, at T = 0, f/^ ~ t/c — 
upon requiring fl ^ Ub one immediately sees that, due to the localized modes of the hidden 
spectrum, the chemical potential on the comb's backbone is smaller than the one measured 
on the chain. 

By substituting the wavefunctions of the eigenstates of the hidden spectrum 
Eqs. ( ITOj) and (fTTj) and using jl ^ Ub one gets 

AbV r'^ „ cos A; , //? \ 

Ab = Ac + ^- / dk ■ tanh eA. (53) 

TT Jq Gfe Vl + cos2 k \^ J 

where 6^= a/A^ + AE"^ (1 + cos^ k). The hidden spectrum eigenstates contribute also to 
the gap function A^ through the second term in the rhs of Eq. (153|) : without an hidden 
spectrum. A;, equals Ac. 
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When > Afe, A^, Eq. ([53]), at T = 0, yields 

Afe(T = 0) 1 



= /C (54) 



Ae(T = 0) 1 _ JZ^ 

where rj^ = (l/>/2) log (l + ^/2). Furthermore, at low temperatures, 

A,(T)/A,(T) ^ A5(T = 0)/A,(T = 0). (55) 

Using the parameters Ej and V obtained from the measurements carried on the JJ chain 
(see Appendix A), one gets /C ~ 1.13. 

Upon requiring that the T = backbone's gap function has a BCS hke functional form- 
i.e., Afc(T = 0) = SEje^'^'^^-'^^ , with Ej and V being the renormalized Josephson energy 
and the renormalized interaction term- one is able to estimate the renormalization of the 
Josephson coupling within the LBdG approach. Namely, one has 

Ej = K,Ej- V = /CV, (56) 

which embodies the effects of the hidden spectrum on the Josephson current. In Fig. 5 



we plot, as a function of the normalized temperature, the values of Ic measured in |21 | 
(squares) and the values of Ic obtained from the Ambegaokar-Baratoff [37] formula using 
both the renormalized coupling given by Eq. fl56l) and the gap function along the backbone 



for the comb-like JJN studied in 



21 



22[ | (solid curve): the agreement of the results of the 



LBDG analysis with the outcomes of experiments is very good at low temperature. 



B. Spatial BEC of bosons hopping on a comb 

The thermodynamic properties of non-interacting bosons hopping on a comb evidences 
the emer gen ce of a spatial BEC on the backbone even if the euclidean dimension of a comb 



is 1 



30 



3l| . In this section we shall review the argument showing how this emergent 
behavior follows from the existence of an hidden part in the spectrum of the adjacency 
matrix describing the comb's connectivity; in fact, the presence of a dense set of states 
filling the gap between —2t and Eq induces a change in the chemical potential, which, in 
turn, allows for the existence of a finite value of the critical temperature at which spatial 
BEC occurs. 
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FIG. 5: Critical currents (in units of the critical current on the reference chain at T = 1.2K) as a 
function of T/Tc for the backbone and the chain. The solid lines are the estimated critical currents 
for the backbone (top) and the chain (bottom). Circles (squares): experimental values for the 
chain (backbone). 



Fixing the number of particles in the grand canonical partition function amounts to 
choose the fugacity z so that 

In Eq. ( !57l) Nt is the number of particles, d{E) is the degeneracy of each single-particle 
eigenstate, /3 = l/fc^T and Eq is the energy of the lowest energy state allowed to a particle 
hopping on a comb; the sum is taken over the entire spectrum a. 

For free bosons hopping on a square comb of size L (with Ns = L x L sites), one has 

Nr = N,,mT) + N.AL,T) + N.,mT) + J^^^^ dE (58) 

where Neo{L,T), N„_{L,T) and N„_^_{L,T) denote, respectively, the number of particles 
which, at a certain temperature T, occupy the ground-state and the two regions cr_ and 0"+ 
of the hidden spectrum; p{E), with E G ctq, is the energy density of states defined in Eq. 
(HOj) . Upon defining tieq = Neq/L"^, ti^q = N^^/L"^ and = N^^/L"^ as the contribution 
given by each one of the three spectral regions to the number of particles per site, one has 

20 



as L — s> oo - that 

nEo{T) = lim / , 



2 ^^"^^^^ 1 



L^oo ^- 1 ^/3t [ v/l+cos2 (27m/L)] 



and 



{L-l)/4 



L^oo ^_ 1 [V8+2^/l+cos2 (27m/L+7r/2)] 1 

n=l ^ e V X 

2 1 

< lim TT^- = V T. 

L^oo L 2-le/3t(y8+2) _ I 

The last equation shows that, at any finite temperature T, ct+ is not macroscopically occu- 
pied. 

The last term of the right-hand side of Eq. fl5S]) is the number of bosons in the delocalized 
states. On a chain the integral appearing in Eq. fl58l) is diverging since the limit z — > 1 is 
attained when the energy equals —2t; at variance, for a comb, the existence of the hidden 
spectrum renders the same integral convergent since the limit z ^ 1 is now attained at the 
lower energy —2t\/2 which lies outside the interval —2t, 2t. One then sees explicitly how the 
network's inhomogeneity works to induce the emergence of a spatial BEG even if the bosons 
are "free" and the network is one-dimensional. 

If one defines as Tc the critical temperature at which spatial BEG on the comb's backbone 
occurs, then, for any T < Tc, the localized ground-state is macroscopically filled. Since 
^EoiTc) = 0, using Eqs. PUI) and (!58l) . Tc may be determined (as a function of the filling 
fraction / and of the hopping strength t) from 

f^' dE 1 

Equation (|59ll can be solved numerically for any value of /. When / ^ 1, one may expand 
the exponential in Eq. fl5^ to the first order in the inverse of the critical temperature 
getting 

- ' dE^====—-r. (60) 



ksTc J.2t V4t2 -E^E-E^ 
Substituting cos6' = El2t in Eq. ( l60l) . one has 

2tTif r de 



ksTc Jo cose-Eo/2t'' 
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from which 



ksTc = Ej 



\ 



2 



2f C'l) 

with Ej being the Josephson energy defined in Eq. fl27|) . Equation (|6T1) allows for a simple 
estimate of the critical temperature Tc at which BEC occurs for free bosons hopping on a 
comb. Upon inserting the pertinent value of the ground-state energy Eq in Eq. fl60|) . one 
finally gets 

ksT, ^ Ej. (62) 

A similar analysis can be carried out for the variety of comb-shaped networks depicted in 
fig. 4. 

The condensate's fraction may be easily determined as a function of the scaled tempera- 
ture T = ^. Taking into account that 

- 2£ 



expanding to the first order in (3 the exponential appearing in Eq. |38| . and, finally, 
using Eqs. (158|) and (!63l) . one may easily compute the number of particles occupying the 
states belonging to a_. Namely, one has to compute A^o = + N^., where Ne^ (^o-_) is 
the number of particles in the ground state (in (T_, except for those occupying the ground 
state): as a result the fraction of condensate, for T < Tc, is given by 

For / ranging from 10^ to 10^, the results provided by Eq. flM|) differ from those obtained 
by the numerical evaluation of Nq from Eq. (I58l) by less than 1%. Equation shows that 
the condensate has dimension 1. 

Finally, one may compute the average number of bosons Nsix^y) occupying a generic 



site (x,y) of a comb 321]; of course, this number depends only on y due to the translational 



invariance of a comb along the backbone. One finds that Nsiy, t) is given by: 

E,.ea^ (65) 

In Eq. (l65l) ipEoiv) is the wavefunction corresponding to the ground-state of the single- 
particle spectrum and tPeSv) the eigenfunctions corresponding to the energies E^ of the 
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hidden spectrum o"_; N„_{En) is the number of particles with energies En G o"- and A^^;^ 
is the number of particles in the ground-state. In the last term of Eq. (l65l) ipEiu) are the 
delocalized eigenf unctions of the eigenvalue equatio n f[29l) . For determining NB{y]T), one 
needs to compute A'"^;;, and Na-_, which are evaluated [32] in Appendix B. Using these results, 
an explicit analytical form for the number of bosons at site y, Nsiy), may, then, be derived. 
The last term in Eq. fl65l) yields, in fact, the contribution coming from the delocalized states: 
for a large network (L ^ 1), and far away from the backbone, this number is independent 
on the site index y and equals a constant {Nt/L'^)t. Using Egs . (|B7D and (IBSp . for r < 1, 
one finds that the distribution of bosons is, for ?/ ^ 1, given by 32] 

„ r. (66) 

The signature of the emerging spatial BEC in a system of non-interacting bosons hopping 
on a comb-shaped network is provided then by the remarkably sharp decrease of the number 
of bosons at sites located away from the backbone. The linear dependence exhibited by the 
solid line in Fig. [6] is consistent with the observation that the condensate has dimension 1. 
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V. CONCLUDING REMARKS 



We evidenced how the optimal engineering of the shape of a network leads to the emer- 
gence of complex features in quantum devices realized with bosonic and superconducting 
networks, whose connectivity is described by an adjacency matrix supporting an hidden 
spectrum. For this purpose we analyzed the paradigmatic case of comb-shaped bosonic and 
superconducting networks. 

For free bosons hopping on a comb, we evidenced how the network's connectivity is 
responsible for the emergence of a spatial BEC along the comb's backbone and computed 
the critical temperature where a spatial BEC emerges. We then analyzed the inhomogeneous 
distribution of the bosons along the comb fingers and estimated the dependence of the non- 
condensate fraction on the reduced temperature r; we showed that the signature for the 
emergence of a spatial BEC on the comb's backbone is provided by a rather sharp decrease 
of the number of bosons occupying the fingers as the temperature is lowered below T^. Finite 
size corrections to our results are negligible already for / ~ 100. With little modifications 



our analysis could be carried out also for diverse network's topologies supporting BEC |31 |. 

We analyzed also superconducting JJNs fabricated on a comb-shaped insulating support. 
We showed that a non perturbative (i.e., induced by the states of the hidden spectrum) 
renormalization of some (i.e., the ones located on the backbone) of the Josephson couplings 
of a comb-shaped JJN is responsible for the observed enhancement of Ic of the Josephson 
junctions located along the comb's backbone. We used an effective theory based on the 
BdG equations since it allows for a simple and rather intuitive derivation of Eq. (1531) , which 
evidences the crucial role played by the hidden spectrum in determining the enhancement 
of the Josephson current along the comb's backbone. The BDG approach relies on a few 
key assumptions; namely, that the eigenfunctions of the BdG equations may be written in 
a tight binding form and that only the electrons close to the Fermi surface contribute to 



determine Ej; once these assumptions are made, one is able to derive Eqs. (iTOl) - ( iTTj) and to 



account for all the dependence on the e 
Ej and V, which may be determined 



the parameters 



2l| on the linear 



ectronic states into the definition o:: 
22l | from the measurements carried in 
chain. Our approach yields a value of the renormalized Josephson coupling of the junctions 
located on the comb's backbone in excellent agreement with the experimental results (see 
fig. 5). Similar phenomena happen for the class [31] of JJNs fabricated on graphs whose 
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adjacency matrix supports an hidden spectrum. 

An alternative way to look at comb-shaped networks is to regard them as a linear chain 
immersed in an environment mimicked by the addition of the fingers. This situation may 



be analyzed using either the Caldeira-Leggett [39|] or - for superconducting devices- the 
electromagnetic environment 40|] models. For Josephson devices this point of view was 
advocated long ago in 41]]. For these devices one expects that the nominal value of the 
Josephson energy Ej of the junctions in the array gets renormalized by the interaction 
with the environment. However, one usually assumes that the effective boundary conditions 
for the quantum fiuctuations of the environment modes do not depend on the Josephson 
couplings or on the network's topology: while this assumption may be perfectly legitimate for 
weak environmental fiuctuations, better care should be used if these fiuctuations are strong 
as it may well happen for one dimensional JJNs. A simple paradigmatic example of a non 
perturbative renormalization of Josephson couplings is given by the simple inhomogeneous 



one-dimensional array analyzed in 
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431], where the source of inhomogeneity is given by 



putting on a site of the linear chain a test junction with a different nominal value of the 
Josephson coupling Ej. Our analysis shows that, for a comb-shaped JJN, the Josephson 
couplings on the backbone get renormalized and that this renormalization is non perturbative 
since the peculiar connectivity of a comb modifies the spectrum of quantum modes living 
on linear chains by the (obviously non-perturbative) addition of an infinite set of localized 
states below the continuum threshold: adding the fingers to a backbone chain is, in fact, 
a topological operation since it amounts to a non trivial change of boundary conditions for 
the Josephson linear chain. It would be interesting to investigate in this perspective also 
systems realized with cold atoms living on pertinent inhomogeneous optical lattices. 

Our analysis provides experimentally testable examples of how the space connectivity 
affects coherent behaviors of physical systems. 
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APPENDIX A: SOLUTION OF LBDG FOR A JJ-CHAIN 

Let us consider a Josephson chain with A^^ grains: i = 1, ■ ■ ■ , Ns and a ^ k, with 
k = 27171/ Ns and n = —Ns/2, ■ ■ ■ , Ns/2 — 1 (periodic boundary conditions are used). Since 
Aij = Sij+i + Sij^i, one gets that the eigenvalues of the matrix eij are 

Ek = -2t cos k + Uc- fi, (Al) 

where U{i) = Uc does not depend on the grain position. Similarly, one sets A{i) = Ac. 
One finds then Uk{j) = Ng^^'^UkC^^^ and Vk{j) = N^^^'^Vke^'^^ , which, in turn, lead to f/| = 
(1/2) [1 + Ek/ Gfc] and if = (1/2) [l _ Ek/ Gfc], where Gfe= ^/A? + From one gets 
1 = (y/2Ns) Y.k tanh (/5 Ek /2) and, for Ns > 1, one finally obtains 

^ = L r -^=^ — t^^h f f + {E-fi + Uc)A . ( A2) 

A BCS-like behavior is obtained with Uc — fi ^ 0: Eq. (1A3P for Ac reads then 

1 = , -^ ^ tanh ( ^V^+E^] . (A3) 

At first sight equation flASp appears quite different from the corresponding BCS gap equation 
([9]), basically due to the factor (1 — ii^^/4t^)~^/^ appearing in the integrand and originating 
from the density of states on the chain. Nevertheless, the behavior is basically BCS. In fact, 
at T = 0, one gets 



where K{x) is the complete elliptic integral of first kind 45|. If one takes the limit Ac/t <C 1 



4*2 



(which is, in a sense, the equivalent of the limit n{0)VBcs 1 in the BCS theory) and recalls 



that, at T = 0, K{x) ^ (1/2) log 16/(1 - x) [isj, one obtains 

Ac(T = 0) ^ 8te-2-*/v^ (A5) 

which is the ID equivalent on a chain of the well-known BCS expression A(T = 0) = 

2;io^^e-i/"(°)^sc's. At T = Tc {Ac{T = Tc) = 0), one has 

V de , / et \ , 

tanh — — , (A6) 



271^^0 eVl^ \2kBT, 
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FIG. 7: A/A(T = 0) vs. T/Tc obtained from the numerical solution of the BCS gap equation ([9]) 
with the Nb parameters (solid line) and from the numerical solution of the gap equation for the ID 
chain of superconducting grains ()A3p . with the parameters obtained for the setup of |21j] (dotted 
line) . 



which, for Ac/t ^ 1, yields 1 ~ {V /2-Kt){\ogt/kBTc + \ogC) with C ~ 4.536. As a result one 
gets 

ksT, = (A7) 

which is the ID equivalent on a chain of the well-known BCS expression ksTc = 
1.14/icjDe"^/"(°)^^^^. Combining Eqs. ^EB and (lA7| enables to show that 

ksTc C 

which coincides with the result ex pec ted from the BCS theory. 

For the experimental setup of 2l|, |22l] one can estimate t ^ ks ■ 70K and V/t ~ 1.75, 
yielding ^ 8.6K and Ac(T = 0) ^ fcs ■ 15.9K. In Fig. Hi we plot the ratio A(r) /A{T = 0) 
vs. the reduced temperature T/T^ obtained from the numerical solution of the BCS gap 
equation ([9]) with the Nb parameters for the bulk (n(O)VBCs = 0.28 and Hud = ks ■ 275K) 
and the one obtained from the numerical solution of the gap equation for the ID chain 
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of superconducting grains (lA3p . with the parameters obtained for the setup of 2l|. The 
comparison evidences the close similarity between the two curves. 

Measurements on a chain made with Nh grains yield ~ and AJT = 0) ~ 
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22j it is 



lAmeV ~ ks ■ 15.9K; furthermore, in the experimental setup described in 
Ic ~ 18^A. The parameters Ej and V, determined from Eq. (lA7p are then given by Ej 
kB ■ ^30K and V/Ej = 1.185. 



APPENDIX B: DISTRIBUTION OF BOSONS IN THE HIDDEN SPECTRUM 



In the thermodynamic limit, for r < 1, Nrj_{En, r) is given by 44l | 

2 



NUE^;r)=lun^L^^ 



L2 



(Bl) 



and it depends on the number of particles in the ground-state Neq- Since fc^Tc ~ 2tf, from 
Eq. (IBip . it follows that the number of particles occupying the hidden spectrum is given by 



(^) = Yl (^«) = -^^0 + 



r Nt 



n=l 



coth 



Eo 



T Nt 



V2N, 



Eo 



(B2) 



As a consequence of the fact that A^o = + Neo, from Eq. (Ell), one finds the way to 
determine Neq/Nt as a function only of the scaled temperature r: 



T 



Eo' 



coth 



r 



A^T 



Eo 



(B3) 



V^Neo 

Solving Eq. (1B3|) and substituting back the value obtained for Neo in Eq. (IBip allows for an 
exact numerical evaluation of Eq. fl65|l . 

Conventional wisdom supported by numerical evidence suggests that - apart from a small 
range of temperatures near - the largest contribution to Nq comes from Neo- Thus, it is 
physically appealing to assume that A'^;;, is given by 



NEo=NT{l-T)g{T). 



(B4) 



In Eq. (IB4[) . gij) is a function only of the scaled temperature r and parametrizes the 
contributions to A"o coming from the states belonging to the hidden spectrum: when g = l, 
the condensate is in the ground-state, while, for = 0, is in the states of the hidden spectrum. 
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FIG. 8: The function g{T) defined in Eq. ()B4p : the empty circles correspond to the numerical solu- 
tion of the self-consistency equation ()B5P : the solid line corresponds to the approximate expression 
(|B6]). 



Substituting Eq. flB4|) in Eq. fiBll) and requiring A^o = + Neq with A^o given by Eq. fl6ll) . 
leads to a self-consistency equation for g{T): 

^ ^ ^ (B5) 



V^(l-r)(7(r) 



coth 



V2{1 - T)g{T) 



1. 



For r not too close to 1, a rather simple approximate solution of Eq. fIBSP is provided by 

(B6) 



coth 



r 



v^(l 



V2{1 - t) 

The error made in using Eq. (1B6P instead of the exact solution of Eq. (IBSp is within few 
percents: for r < 0.5 the error is less than 1%, while for r = 0.7 is about 5%. In Fig. (jS]) 
we plot the function g{T) as obtained from the numerical solution of the self-consistency 
Eq. (IBSP and from the approximate expression (]B6p . 
Upon inserting Eq. flB6p in Eq. flB4p . one has 




V2(l - r) 



coth 



V2(l 



and, from Eq. flBlP . one gets 
iV,_(r) ^ 2 



E 



V2(7rn)2 
n=l T r 



(1- 




coth 



V2{1 



(B7) 



[B8) 



[1] N. Goldenfeld and L. Kadanoff, Science 284, 87 (1999). 



29 



[2] E. Dagotto, Science 309, 257 (2005). 

[3] P. Xiong et al., Phys. Rev. Lett. 69, 3220 (1992); E. Fradkin and S. A. Kivelson, Phys. Rev. 

59, 8065 (1999); S. A. Kivelson, E. Fradkin and V.J. Emery, Nature 393, 550 (1998). 
[4] E. Dagotto, T. Hotta and A. Moreo, Phys. Rep. 344, 1 (2001); E. Dagotto, Nanoscale Phase 

Separation and Colossal Magnetoresistance, Springer- Verlag (2001). 
[5] V.J. Emery, S.A. Kivelson and J.M. Tranquada, Proc. Natl. Acad. Sci. U.S.A. 96, 8814 (1999); 

J. Zaanen, Nature 404, 714 (2000) and references therein. 
[6] Ch. Renner et al., Phys. Rev. Lett. 80, 149 (1998); M. Suzuki, T. Watanabe and A. Matsuda, 

Phys. Rev. Lett. 82, 5361 (1999); T. Ekino et al. J. Low Temp. Phys. 117, 359 (1999). 
[7] S.A. Kivelson and E. Fradkin, Preprm^ |cond-mat/0507459 

[8] Yu.N. Ovchinnikov, S.A. Wolf and V.Z. Kresin, Phys. Rev. B63, 064524 (2001). 

[9] P.W. Anderson, Science 177, 393 (1972). 
[10] see, for instance: Coherence in Superconducting Networks, J.E. Mooij and G. Schon eds., 
Physica B152, pp. 1-308, (1988); Josephson Junction Arrays, H.A. Cerdeira and S.R. Shenoy 
eds., Physica B222, pp. 253-406, (1996); R. Fazio and H. van der Zant, Phys. Rep. 355, 235 
(2001). 

[11] O. Morsch and M. K. Oberthaler, Rev. Mod. Phys. 78, 179 (2006); F.S. Cataliotti et al, 

Science 293, 843 (2001); B.P. Anderson and M. Kasevich, Science, 282, 1686 (1998). 
[12] E. Simanek, Inhomogeneous Superconductors, Oxford University Press, New York, 1994. 
[13] Y. Makhlin, G. Shoen and A. Shnirman, Rev. Mod. Phys. 73, 357 (2001). 
[14] P.G. de Gennes, C. R. Acad. Sci. Ser. B292, 279 (1981); P.G. de Gennes, C. R. Acad. Sci. 

Ser. B292, 9 (1981); S. Alexander, Phys. Rev. B27, 1541 (1983); H.J. Fink, A. Lopez, and R. 

Maynard, Phys. Rev. B26, 5237 (1982); R. Rammal, T.C. Lubensky, and G. Toulouse, Phys. 

Rev. B27, 2820 (1983). 

[15] G. Deutscher and R. Rosembaum, Appl. Phys. Lett 27, 366 (1975); G. Deutsher, L Grave and 
S. Alexander, Phys. Rev. Lett. 48, 1497 (1982); G. Deutscher et al, Phys. Rev. B24, 6464 
(1981). 

[16] J. Berger, J. Rubinstein eds. Connectivity and Superconductivity, Lecture Notes in Physics, 

Springer- Verlag, Berlin, 2000. 
[17] L.B. loffe et al. Nature 415, 503 (2002); B. Doucot, M.V. Feigel'man and L.B. loffe, Phys. 

Rev. Lett. 90, 107003, (2003); B. Doucot, L.B. loffe and J. Vidal, Phys. Rev. B69, 107003 

30 



(2003); B. Doucot et al. Phys. Rev. B 71, 024505 (2005). 
[18] X.G. Wen and Q. Niu, Phys. Rev. B41, 9377 (1990); X.G. Wen, Phys. Rev. Lett. 90, 016803 

(2003) . 

[19] C.C. Abiho et al, Phys. Rev. Lett. 83, 5102 (1999); J. Vidal, R. Mosseri and B. Doucot, Phys. 

Rev. Lett. 81, 5888 (1998). 
[20] I.V. Protopopov and M. V. Feigel'man, Phys.Rev. B70, 184519 (2004); I.V. Protopopov and 

M.V. Feigel'man, |cond-mat/0510766[ 
[21] P. Silvestrini, R. Russo, V. Corato, B. Ruggiero, C. Granata, S. Rombetto, M. Russo, M. 

Cirillo, A. Trombettoni, and P. Sodano, cond-mat/05 12478 
[22] P. Sodano, A. Trombettoni, P. Silvestrini, R. Russo, and B. Ruggiero, New J. Phys., 8, 327 

(2006). 

[23] L. P. Pitaevskii and S. Stringari, Bose-Einstein Condensation, Oxford University Press, Ox- 
ford, 2003. 

[24] M. K. Oberthaler and T. Pfau, J. Phys.: Condens. Matter 15, R233 (2003). 
[25] R. Roth and K. Burnett, Phys. Rev. A 68, 023604 (2003). 

[26] I. Brunelli, G. Giusiano, F. P. Mancini, P. Sodano, and A. Trombettoni, J. Phys. B 37, S275 

(2004) . 

[27] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler, Phys. Rev. 

Lett. 95, 010402 (2005). 
[28] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Phys. Rev. Lett. 79, 4950 (1997). 
[29] F. Harary, Graph Theory, Addison- Wesley, Reading (1969). 

[30] R. Burioni, D. Cassi, I. Meccoli, M. Rasetti, S. Regina, P. Sodano, and A. Vezzani, Europhys. 
Lett. 52, 251 (2000). 

[31] R. Burioni, D. Cassi, M. Rasetti, P. Sodano, and A. Vezzani, J. Phys. B 34, 4697 (2001). 
[32] G. Giusiano, F. P. Mancini, P. Sodano, and A. Trombettoni, Int. J. Mod. Phys. B 18, 691 
(2004). 

[33] P.G. de Gennes, Superconductivity of Metals and Alloys, Addison- Wesley (1989). 
[34] B. Abeles, Phys. Rev. B 15, 2828 (1997). 

[35] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 
(1998). 

[36] F. P. Mancini, P. Sodano, and A. Trombettoni, to appear in "Eleventh Training Course in 

31 



the Physics of Correlated Electron Systems and High-Tc Superconductors", Vietri sul Mare, 

Italy, Oct 2006, |cond-m at/061238 8[ 
[37] V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 (1963); ibid. 11, 104 (1963) 
[38] This approximation holds for / S> 1 and it is in very good agreement with the results of a 

numerical evaluation of the integral (j63p . 
[39] A.O. Caldeira and A.J. Leggett, Ann.Phys. (N.Y.) 149, 374 (1983). 

[40] M.H. Devoret et al, Phys. Rev.Lett. 64, 1824 (1990); S.M. Girvin et al, 64, 3183 (1990); G. 

Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 
[41] A. Schmid, J. Low Temp. Phys. 49, 609 (1982). 
[42] L.L Glazman and A.L Larkin, Phys. Rev. Lett. 79, 3736 (1997). 
[43] D. Giuhano and P. Sodano, Nucl. Phys. B711, 480 (2005). 

[44] P. Buonsante, R. Burioni, D. Cassi, and A. Vezzani, Phys. Rev. B 66, 094207 (2002). 
[45] M. Abramowitz and L A. Stegun, Handbook of mathematical functions with formulas, graphs, 
and mathematical tables, National Bureau of Standards, Washington, D.C., 1964. 



32 



